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We consider the possibility of a gravitationally induced particle production through the mechanism 
of a nonminimal curvature-matter coupling. An interesting feature of this gravitational theory is 
that the divergence of the energy-momentum tensor is nonzero. As a first step in our study we 
reformulate the model in terms of an equivalent scalar-tensor theory, with two arbitrary potentials. 

By using the formalism of open thermodynamic systems, we interpret the energy balance equations 
in this gravitational theory from a thermodynamic point of view, as describing irreversible matter 
creation processes. The particle number creation rates, the creation pressure, and the entropy 
production rates are explicitly obtained as functions of the scalar field and its potentials, as well as 
of the matter Lagrangian. The temperature evolution laws of the newly created particles are also 
obtained. The cosmological implications of the model are briefly investigated, and it is shown that 
the late-time cosmic acceleration may be due to particle creation processes. Furthermore, it is also 
shown that due to the curvature-matter coupling, during the cosmological evolution a large amount 
of comoving entropy is also produced. 

PACS numbers: 04.50.Kd,04.20.Cv, 95.35.+d 


I. INTRODUCTION 


The simplest explanation of the late-time accelerated 
expansion of the Universe [ 2 , [ 2 } is to invoke a cosmological 
constant, A, which can be associated to the vacuum en¬ 
ergy i- Despite the excellent fit to observational data, 
the presence of A suffers from two serious drawbacks, 
namely, the cosmological constant problem and the coin¬ 
cidence problem. However, it has recently been argued 
that the vacuum energy is not constant but decays into 
other particle constituents [ij. Indeed, phenomenolog¬ 
ical models, with a variable cosmological constant [5|, 
have been proposed to address the above problems. For 
instance, a simple and thermodynamically consistent cos¬ 
mology with a phenomenological model of quantum cre¬ 
ation of radiation due to vacuum decay was presented 
in [6j, where the thermodynamics and Einstein’s equa¬ 
tions lead to an equation in which H is determined by 
the particle number N. It was shown that the evolution 
equation for H has a remarkably simple exact solution, 
in which a non-adiabatic inflationary era exits smoothly 
to the radiation era, without a reheating transition. In 
[7], a new accelerating flat model without dark energy 
that is fully dominated by cold dark matter (CDM) was 
investigated. It was shown that the number of CDM par- 
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tides is not conserved and the present accelerating stage 
is a consequence of the negative pressure describing the 
irreversible process of gravitational particle creation. 

In [f[ , the correspondence between cosmological mod¬ 
els powered by a decaying vacuum energy density and 
gravitationally induced particle production was explored. 
Although being physically different, it was shown that 
under certain conditions both classes of cosmologies can 
exhibit the same dynamical and thermodynamical be¬ 
havior. By using current type la supernovae data, recent 
estimates of the cosmic microwave background shift pa¬ 
rameter and baryon acoustic oscillations measurements, 
the authors performed a statistical analysis to test the 
observational viability of the models and the best-fit of 
the free parameters was also obtained. Furthermore, the 
particle production cosmologies (and the associated de¬ 
caying A(t)-models) were modelled in the framework of 
field theory by a phenomenological scalar field model. 

In this context, a new cosmic scenario wit h g ravitation- 
ally induced particle creation was proposed [8[ , where the 
Universe evolves from an early to a late time de Sitter 
era, with the recent accelerating phase driven only by 
the negative creation pressure associated with the cold 
dark matter component. The model can be interpreted 
as an attempt to reduce the so-called cosmic sector (dark 
matter plus dark energy) and relate the two cosmic ac¬ 
celerating phases (early and late time de Sitter expan¬ 
sions) . A detailed thermodynamic analysis including pos¬ 
sible quantum corrections was also carried out. For a 
very wide range of the free parameters, it was found that 
the model presents the expected behavior of an ordinary 
macroscopic system in the sense that it approaches ther- 
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modynamic equilibrium in the long run (i.e., as it nears 
the second de Sitter phase). Moreover, an upper bound 
was found for the Gibbons-Hawking temperature of the 
primordial de Sitter phase 9]. Finally, when confronted 
with the recent observational data, the current ‘quasi’-de 
Sitter era, as predicted by the model, it was verified to 
pass the cosmic background tests very comfortably. 

In this work, we consider an alternative mechanism for 
the gravitational particle production, namely, through a 
nonminimal curvature-matter coupling, in modified the¬ 
ories of gravity. A general property of these theories 
00 ], is the non-conservation of the energy-momentum 
tensor (for a recent review of modified gravity models 
with curvature-matter coupling see 0 )- Thus, the cou¬ 
pling between the matter and the higher derivative cur¬ 
vature terms may be interpreted as an exchange of energy 
and momentum between both. This latter mechanism in¬ 
duces a gravitational particle production. We note that 
the generalized energy balance equations in these grav¬ 
itational theories have been interpreted from a thermo¬ 
dynamic point of view as describing irreversible matter 
creation processes in 0- Thus, the coupling between 
matter and geometry generates an irreversible energy 
flow from the gravitational field to newly created mat¬ 
ter constituents, with the second law of thermodynamics 
requiring that the geometric curvature transforms into 
matter. 

Here we extend and refine the analysis initiated in 0 
by investigating in detail the thermodynamic interpre¬ 
tation of the curvature-matter gravitational coupling for 
the so-called linear version of the f(R,L m ) gravity the¬ 
ory, where R is the Ricci scalar and L m is the mat¬ 
ter Lagrangian, which we denote as Lf(R , L m ) theory, 
with the gravitational Lagrangian given by fi(R)/2 + 
[1 -I- A /2 (R)} L m , where fi(R) and f 2 {R) are arbitrary 
functions of R , and A is a coupling constant. We note 
that in the Lf(R,L m ) theory, the gravitational action is 
linear in the matter Lagrangian L m , and not in the Ricci 
scalar R. As a first step in our study we introduce the 
equivalent scalar-tensor description of the theory [l8[, in 
which the action is equivalent to a Brans-Dicke type the¬ 
ory, with a single scalar field if), a vanishing Brans-Dicke 
parameter w, and a coupling between the scalar 

field and matter. 

By using the formalism of open thermodynamic sys¬ 
tems Il9l-l22|. we interpret the energy balance equation 
of the theory as describing a matter creation process. In¬ 
deed, the irreversible thermodynamics of open systems, 
and its implications for cosmology have been extensively 
analyzed 1231] . Here, we obtain the equivalent particle 
number creation rates, the creation pressure and the en¬ 
tropy production rates as functions of the scalar field, of 
the two scalar potentials, and of the matter Lagrangian, 
respectively. The temperature evolution of the newly cre¬ 
ated particles is also obtained. Due to the curvature- 
matter coupling, during the cosmological evolution a 
large amount of comoving entropy could be produced. 
The cosmological implications of the theory in its scalar- 


tensor representation are also investigated. 

The present paper is organized as follows. In Section 
im we present the action and the field equations of the 
modified gravity model with a linear curvature-matter 
coupling in both their standard and scalar-tensor rep¬ 
resentations. The thermodynamic interpretation of the 
theory is developed in Section Hill where the particle cre¬ 
ation rates, the creation pressure and the entropy pro¬ 
duction are analyzed in detail. In Section IIV1 the cos¬ 
mological implications of the theory in its scalar-tensor 
representation are considered. In SectionfVl the behavior 
of the entropy of the Universe in the Lf ( R , L m ) gravity 
theory, with the horizon entropy included, is analyzed. 
We discuss and conclude our results in Section IVTl 

II. NONMINIMAL CURVATURE-MATTER 
COUPLING 

A. General formalism 

The action of f(R) gravity can be generalized with 
the introduction of a linear nonminimal coupling between 
matter and curvature. The corresponding gravitational 
theory, a particular case of the general f(R,L m ) theory 
0 , and which we denote as Lf ( R , L m ), has the action 
given by 0 , 

S = j | + [1 + -VM-R)] L m (.9fivi $) j V~9 d A x 

(1) 

where the factors fi(R ) (with i = 1,2) are arbitrary func¬ 
tions of the Ricci scalar R. The coupling constant A 
determines the strength of the interaction between f 2 (R) 
and the matter Lagrangian. L m is the matter Lagrangian 
density, which is a function of the metric and of the 
matter fields $. 

Now, varying the action with respect to the metric g 
provides the following field equations: 

QRiiv - \fi(R)g^ + = [1 + A f 2 {R)]T liV , ( 2 ) 

where F^R) = f-{R), i = 1,2, and the prime represents 
a derivative with respect to the scalar curvature R. We 
have defined 


Q = F 1 {R) + 2XF 2 (R)L m , 

( 3 ) 

= (g^O - V M V„), 

( 4 ) 


for notational simplicity. The matter energy-momentum 
tensor is defined as 

rp _ _ 2 5(^jL m ) 

V=9 S(g^) ■ [ j 

An important property of any gravitational theory 
is its stability with respect to local perturbations. 
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In the standard f(R) gravity, a fatal instability (the 
Dolgov-Kawasaki instability) appears once the condition 
f"(R)< 0 0. The instability develops on time scales of 
the order of 10 26 s. The stability properties of the gravi¬ 
tational models described by the action given by Eq. 0 
were studied in j25|, and it turns out that the correspond¬ 
ing stability criterion is (R) + 2A f 2 {R) > 0. In order 
to obtain the stability condition one expands the param¬ 
eters of the model as the sum of a background field with 
constant curvature, and a small perturbation, so that 
R = Ro + Rh T = T 0 + T 1 , fi (R) = Rq + I?i + ecj) (Ro) + 
ecj)' (Rq) R\ + ..., f[(R) = 1 + ccj)' (Rq) + e4>" ( Rq ) Ri + .... 
With the use of the linearized field equations we obtain 
the stability condition e(j>" (R) + 2Xf 2 (R) > 0 [2f|, gener¬ 
alizing the stability!condition f"(R) = ecj>"(R) > 0, found 
in f(R) gravity [24j . From a physical point of view the 
reason for stability is once the stability conditions are 
satisfied the effective mass m e // of the dynamical degree 
of freedom associated to the small perturbation R\ of the 
background curvature is non-negative. 

A general property of these nonmininral curvature- 
matter coupling theories is the non-conservation of the 
energy-momentum tensor. This can be easily verified by 
taking into account the covariant derivative of the field 
Eq. @, the Bianchi identities, = 0, and the fol¬ 

lowing identity, 

(□V, - V„D)Fi = R„ v WFi, (6) 


An intriguing feature is that the extra force depends 
on the form of the Lagrangian density. Note that con¬ 
sidering the Lagrangian density L m = —p, where p is the 
pressure, the contribution of the nonminimal curvature- 
matter vanishes [28] . It has been argued that this is not 
the unique choice for the matter Lagrangian density and 
that more natural forms for L m , such as L m = p, do not 
imply the vanishing of the extra-force. Indeed, in the 
presence of the nonminimal coupling, they give rise to 
two distinct theories with different predictions 271*30. 

B. Scalar-tensor representation of the linear 
curvature-matter coupling 

As has been shown in f(R) gravity is equiva¬ 

lent to a scalar-tensor theory. In this context, the equiv¬ 
alence between the modified gravity models to a lin¬ 
ear curvature-matter coupling and scalar-tensor gravity 
models was also established in [T||]. More specifically, it 
was shown that the action given by Eq. m is equivalent 
with a two-potential scalar-tensor Brans-Dicke type the¬ 
ory, with a single scalar field, a vanishing Brans-Dicke 
parameter ui, and an unusual coupling of the second po¬ 
tential U(i/j) of the theory to matter. 

As a first step in the scalar-tensor formulation of the 
theory we introduce a new field 0, and reformulate the 
action m as 


which then implies the following relationship: 

XFo 

V = — -f- [g^L m - T^] W^R. (7) 
I + A/2 

Note that in the absence of the coupling, A = 0, one ob¬ 
tains the conservation of the energy-momentum tensor 
[27| . which can also be verified through the diffeomor- 
phism invariance of the matter part of the action. The 
conservation of the energy-momentum tensor also follows 
from Eq. 0, if f 2 (R) is a constant or the matter La¬ 
grangian is not an explicit function of the metric. 

In order to test the motion in our model, we consider 
for the energy-momentum tensor of matter a perfect fluid 


— (p F p) UfJJu PQfivj (8) 


where p is the total energy density and p, the pressure, 
respectively. The four-velocity, C/ M , satisfies the condi¬ 
tions UfJJ 11 = 1 and = 0. We also introduce 

the projection operator h^x = g^x — U^Ux from which 
one obtains h^xU^ = 0. 

From Eq. 0, we deduce the equation of motion for a 
fluid element: 


DU a 

ds 


dU a 

ds 


+ r °u*u v 


= r 


where the extra force is given by: 


( 9 ) 


S — J ^ 2^ + 2 (-^ — <^)+[l + ^h( 4 >)} Am|- 

( 11 ) 

Next, we introduce the second field ^((j)) = (with a 

prime denoting a differentiation with respect to ^>), and 
thus we obtain for the action the final expression 



V(ip) +U(tp)L m , 


( 12 ) 


where the two potentials V(tp) and U(il>) of the theory 
are defined as 


= <t>(^)fi [</>M] ~ h [</>W] 


(13) 


and 


UW = 1 + A/a [<t>m , (14) 

respectively. The function <fi(il>) must be obtained by 
inverting ip((t>) = /((</>). The actions 0 and (fl2l) are 
equivalent when f\(R) A 0 [l8|, similarly to the case of 
pure f(R) gravity |3lH34j . 

The action given by Eq. m can be written, via a con¬ 
formal rescaling gp v = exp (—a<j>/ 2) g^, a = constant, as 
a four-dimensional dilaton gravity whose action, in the 
“Einstein frame”, has the form [25| 


r = 


i 


XFo 


p + p [l + Xf 2 


(Lm + p) V U R + V„p 


h av ( 10 ) 


S F = 


J d 4 xy/=g e . (15) 
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In the Einstein frame representation of the modified 
gravity with linear coupling between matter and geom¬ 
etry the extra force is due to the coupling between the 
matter Lagrangian and the Brans-Dicke-like scalar (p. On 
the other hand it is important to mention that in the 
Lf ( R , L m ) theory there is no scaling of units with some 
powers of the conformal factor of the conformal transfor¬ 
mation [25| . For this reason it is impossible to reduce 
Lf ( R , L m ) gravity to a standard scalar-tensor theory, or 
to find a string gravity equivalent. 

By varying the action m with respect to g^ u provides 
the gravitational field equations of the scalar-tensor the¬ 
ory as 


that the energy-momentum tensor has the perfect fluid 
form given by Eq. ©, then Eq. m can be written in 
the equivalent form 


(V'V + V» UpUv + {p + p) U v WUp - Wpg^ 
+ (p + p) + V" [In U (if)} T„ v 

21/' (VO - U'(iP)L r 


Uty) 


V„^ = 0. (23) 


By multiplying Eq. (l23l) with U v we obtain the energy 
balance equation in the scalar-tensor representation of 
the linear curvature-matter coupling given by 


V 7 ^ ~ ^g^Rj + P^ip = U(ip)T,j V + V(ip)g,j, v , (16) 

while the variation of the action with respect to the field 
ip gives the relation 

§ - E'(V0 + U'(4>)L m = 0. (17) 

The contraction of the field equation Eq. (HU) yields 
the scalar relation 

- ipR + 3[3ip = U(ip)T + W(ip), (18) 


p + 3 H(p + p) + p— In U (ip) + 
as 


2V' (iP) - U'(ip)L m 
U(iP) 


ip = 0, 
(24) 


where we have introduced the Hubble function H = 
(1/3)V M (7 /i , and we have denoted ‘ = = d/ds, 

respectively, where ds is the line element corresponding 
to the metric g^, ds 2 = g^dx^dx 1 '. After acting on 
Eq. (l23l) with the projection operator h uX , provides the 
momentum balance equation for a perfect fluid as 




d 2 x a 

—— + r“ - f a = o 

ds 2 » v J 


(25) 


where T = Tjf is the trace of the energy-momentum ten¬ 
sor. By combining Eqs. 03 and (fl8l) we obtain the field 
equation of the field ip as 

nip = ^u(ip)T+^v(ip)+^ipv'(ip)-^ipu'(ip)L m . (19) 

By eliminating the term Oip, and taking into account 
Eq. d® the field equations (fTfil) provides 


iP 




- = U(ip ) 




where the extra-force is given by 


/“ = h» a 


In U (ip) 


2V' (ip) - U'(ip)L m 

U(ip) 




(26) 

It is important to note that in the scalar-tensor rep¬ 
resentation of modified gravity with a linear curvature- 
matter coupling, the extra-force acting on test fluids is 
non-zero independently of the choice for the matter La¬ 
grangian. 


j V(iP) + \ipV\ip) ~ \pU’(iP)L r 


g^ 


(20) III. GRAVITATIONALLY INDUCED PARTICLE 

CREATION 


Now, taking the covariant divergence of the field equation 
chd, with the use of Eq. ©, we obtain first 


f+ v' (ip) 


\7„ip 


U’ (ip)VnipTjf 
+U (ip) SJnTjf. (21) 


Then, by eliminating R/2 with the help of Eq. (fT71) we 
obtain for the divergence of the energy-momentum tensor 


V pTjf = — [V M In U (V 1 )] Tjf 


2V' (if) - U'(iP)L m 
U(ip) 


Vvlp. 


( 22 ) 

Equation (1221) allows the formulation of the energy and 
momentum balance equations in the scalar-tensor repre¬ 
sentation of the modified theory of gravity with a linear 
coupling between matter and geometry. By assuming 


In the present Section, we analyze the physical inter¬ 
pretation of the curvature-matter coupling in the scalar- 
tensor representation by adopting the point of view of 
the thermodynamics of open systems, in which matter 
creation irreversible processes may take place at a cos¬ 
mological scale MM- A s we have already seen in 
the previous Section, the energy conservation equation 
of the curvature-matter coupling, given by Eq. ([HD , con¬ 
tains, as compared to the standard adiabatic conserva¬ 
tion equation, an extra term, which can be interpreted 
in the framework of the open thermodynamic systems as 
an irreversible matter creation rate. According to irre¬ 
versible thermodynamics, matter creation also represents 
an entropy source, generating an entropy flux, and thus 
leading, in the presence of the curvature-matter coupling, 
to a modification in the temperature evolution. 
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In the following, we investigate only the case in 
which all the non-diagonal components of the energy- 
momentum tensor of the matter are equal to zero, so 
that = 0, ji v. Generally, the energy-momentum 
tensor for a viscous dissipative fluid in the presence of 
heat conduction is given by T = (p + p + II) — (p + 
II)3/jj/ + q^u^ + where II is the bulk viscous 

pressure, q u is the heat flux, and 7r M „ is the tensor of 
viscous dissipation [26| . q^ and must satisfy the con¬ 
ditions q^u 11 = 0 and 7= 7r(( = 0, respectively. In 
the following we neglect the viscous effects in the cosmo¬ 
logical fluid, thus assuming II = 0 and = 0. In a co¬ 
moving reference frame with u M = (1,0), the form of the 
heat flux vector is fixed by the normalization condition as 
q^ = (0 ,q). Therefore in a comoving reference frame all 
components of the form u^q^ of the energy-momentum 
tensor are identically equal to zero, and is a diagonal 
tensor. From the point of view of the thermodynamics of 
the irreversible processes, this condition implies the im¬ 
possibility of heat transfer in the considered gravitational 
system. In particular, this condition is always satisfied 
in homogeneous and isotropic cosmological models, de¬ 
scribed by the Friedmann-Robertson -Walker geometry, 
since in these models the condition To; = 0, i = 1,2,3 
must always hold. 


A. Matter creation rates and the creation pressure 


We assume that the cosmological metric is given by the 
flat isotropic and homogeneous Friedmann-Robertson- 
Walker (FRW) metric, 

ds 2 = dt 2 - a 2 {t) ( dx 2 + dy 2 + dz 2 ) , (27) 


where a(t) is the scale factor, describing the expansion of 
the Universe. In this geometry the cosmological mat¬ 
ter is comoving with the cosmological expansion, and 
therefore the four velocity of the cosmological fluid is 
= (1,0,0,0), while the Hubble function takes the 
form H = d/a, since = '= d/dt. 

To investigate the thermodynamical implications at 
the cosmological scale with a curvature-matter coupling 
we consider that the Universe contains N particles in a 
volume V, with an energy density p and a thermody¬ 
namic pressure p, respectively. For such a cosmological 
system, the second law of thermodynamics, in its most 
general form, is given by [20] 


d 

dt 


(p« 3 )+*4 


dt 


dQ 

dt 


P + P d 
n dt 


(na 3 ) , 


(28) 


where dQ is the heat received by the system during time 
dt , and n = N/V is the particle number density, respec¬ 
tively. Due to our choice of the geometry of the Universe, 
and of the cosmological principle, only adiabatic transfor¬ 
mations, defined by the condition dQ = 0, are possible. 
Therefore in the following we ignore proper heat trans¬ 
fer processes in the Universe. However, as one can see 


from Eq. (128[) . under the assumption of adiabatic trans¬ 
formations the second law of thermodynamics contains 
the term [(p+p)/n]d (na 3 ) /dt 1 which explicitly takes into 
account the time variation of the cosmological particles 
in a given volume V. Hence, in the irreversible thermo¬ 
dynamics description of open systems, even for adiabatic 
transformations dQ = 0, one can consider the “heat” (in¬ 
ternal energy), received/lost by the system, and which is 
due to the change in the particle number n. For modified 
gravity with a curvature-matter coupling the change in 
the particle number is due to the transfer of energy from 
gravity to matter. Thus, via matter creation, gravity acts 
as a source of internal energy, and of entropy. For adi¬ 
abatic transformations dQ/dt = 0, we can reformulate 
Eq. (l28l) in an equivalent form as 

p + 3 (p + p)H = E.—E. (fa q- 3 Hn). (29) 

n 

Therefore, from the point of view of the thermodynam¬ 
ics of open systems, Eq. m, giving the energy balance 
in the presence of a curvature-matter coupling, can be 
interpreted as describing particle creation in an homoge¬ 
neous and isotropic geometry, with the time variation of 
the particle number density obtained as 

h + 3n7J = Fn, (30) 


where the particle creation rate T is a non-negative quan¬ 
tity defined as 


r = - 


i 


p + p I at 


4ln£W) + 


ua >) 


( 31 ) 

Therefore, the energy conservation equation can be re¬ 
formulated in the alternative form 


p + 3(p + p)H = (p + p)T. (32) 

As proven initially in j20j . for adiabatic transforma¬ 
tions Eq. (1281) , describing irreversible particle creation in 
an open thermodynamic system, can be rewritten as an 
effective energy conservation equation, 

^ (pa 3 ) +{p + Pc) ^a 3 = 0, (33) 

or, in an equivalent form, as, 

p + 3(p + p+p c )H = 0, (34) 


where we have introduced a new thermodynamic quan¬ 
tity, p c , denoted the creation pressure and defined as [20] 


Pc 


p+pd (na 3 ) 
n da 3 


P + P 

3 nH 


(h + 3 nH) = 


p + p r 
3 H' 


(35) 


Therefore in modified gravity with a linear curvature- 
matter coupling the creation pressure is given by 


Pc ~~ 


1 

+H 


P dt lnU ^ J ) + 


2V 7 ty) - U'(^)L m 

u(i 0 



( 36 ) 
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Note that from Eq. CD. the coupling between the 
matter and the higher derivative curvature terms may 
be interpreted as an exchange of energy and momentum 
between both. In the standard formulation of the lin¬ 
ear curvature-matter coupling, by taking into account a 
FRW background, and from Eq. CD we obtain the energy 
balance equation as 


p + 3H(p+p) 


A F 2 (R) 

1 + A/ 2 (.R) 


1 ) P R. 


(37) 


Hence, by considering that the mechanism for gravi¬ 
tational particle production is through the nonminimal 
curvature-matter coupling, so that comparing Eqs. ©D 
and (l37)l . we have for the creation pressure 


Pc = 


A F 2 {R) 1 

1 +A f 2 {R) 3 H 


(1 — a)pR , 


(38) 


with the requirement that p c be negative. We will only 
consider the case of L m = —p, i.e., a = —w, in order to 
have a non-vanishing creation pressure. 


B. Entropy and temperature evolution 

According to the basic principles of the thermodynam¬ 
ics of open systems the entropy change consists of two 
components: the entropy flow term d e S , and the entropy 
creation term diS. The total entropy S of an open ther¬ 
modynamic system can be represented as [01 . |20[ 


dS = d e S + diS , (39) 

where by definition diS > 0. Both the entropy flow and 
the entropy production can be obtained from the total 
differential of the entropy given by [20| . 

Td (sa 3 ) = d (pa 3 ) + pda 3 — pd (na 3 ) , (40) 

where T is the temperature of the open thermodynamic 
system, s = S/a 3 is the entropy per unit volume, and p 
is the chemical potential, defined as 

pn = h — Ts, (41) 


where h = p + p is the enthalpy of the system. 

In the case of a closed thermodynamic system and for 
adiabatic transformations we have dS = 0 and diS = 0. 
However, in the presence of a curvature-matter coupling, 
leading to effective matter creation, there is a non-zero 
contribution to the total entropy. For a homogeneous and 
isotropic Universe the entropy flow term d e S vanishes, so 
that d e S = 0. On the other hand matter creation also 
represents a source for entropy creation, and the time 
variation of the corresponding entropy is obtained as [20j 



dS d / q\ d / 

T ~it ( n “ 3 ) - 
n at v 7 


( 42 ) 


Equation ([42]) gives the time variation of the entropy 
as 

^ = - (n + 3 Hn) = TS > 0, (43) 

dt n 

so that the entropy increase due to particle production 
yields the expression 

S(t) = S 0 eJ'o r ( t, ) dt \ (44) 

where Sq = 5(0) is a constant. With the use of Eq. (1431) . 
we obtain for the entropy creation in the scalar-tensor 
representation of the linear coupling between matter and 
geometry the following equation 


ldS_ _ _ 1_ 

S dt p + p 1 r dt 


P Uu W+ 2V '^- u ' WLm A. 


U(i o 


(45) 


The entropy flux four-vector 5 M is defined as |2X] 

5^ = naU 11 , (46) 


where a = S/N is the specific entropy per particle. 5 M 
must satisfy the second law of thermodynamics, which 
imposes the constraint V M 5 M > 0 on its four-dimensional 
divergence. The Gibbs relation [2l| , 

nTda = dp - dn, (47) 

n 

together with the definition of the chemical potential p 
of the open thermodynamic system, 

p=-~Ta, (48) 

n 

yields 

= {h + 3nH)a + nU^\/f,a 
= ^(h + 3 Hn) (^- p'j , (49) 

where we have used the relation 

nTa = p — = 0, (50) 

n 

which immediately follows from Eq. (OHl) . With the use 
of Eq. (l30l) we obtain the entropy production rate due to 
the particle creation processes given by 






T{p + p) 


p^lnUU) 


+ 


2V' (VO - U'W)L m 

u(i 0 




(51) 


A general thermodynamic system is described by two 
fundamental thermodynamic variables, the particle num¬ 
ber density n, and the temperature T, respectively. If 









7 


the system is in an equilibrium state, the energy den¬ 
sity p and the thermodynamic pressure p are obtained, 
in terms of n and T, from the equilibrium equations of 
state of the matter, 


P = p(n,T), p = p(n,T). (52) 


Therefore the energy conservation equation (El can be 
obtained in the following general form 

^h+^f + 3(p + p) H = Tn. (53) 

By using the general thermodynamic relation [2lj] 


dp h T dp 
dn n n dT' 


(54) 


it follows that the temperature evolution of the newly 
created particles due to the curvature-matter coupling is 
given by the expression 


T r> 

= c? s -=cl(T-3H), (55) 

1 n 

where the speed of sound c s is defined as (? s = dp/dp. 
If the geometrically created matter satifies a barotropic 
equation of state of the form p = (7 — 1 ) p, 1 < 7 < 2 , 
the temperature evolution follows the simple equation 

T = Tori' 1 - 1 . (56) 


C. Bulk-viscosity description of matter creation 
processes with a curvature-matter coupling 


An alternative physical interpretation of particle cre¬ 
ation processes in cosmology was suggested by Zeldovich 
[35]. and later on by Murphy [36} and Hu (37}- Accord¬ 
ing to this interpretation, the viscosity of the cosmologi¬ 
cal fluid represents a phenomenological description of the 
effect of the creation of particles by the non-stationary 
gravitational field of the expanding universe. Therefore, 
from a physical point of view, a non-vanishing particle 
production rate is equivalent to the introduction of a bulk 
viscous pressure in the energy-momentum tensor of the 
cosmological fluid. From a quantum mechanical point 
of view, such a viscous pressure can also be related to 
the viscosity of the vacuum [351 - 137} . This physical inter¬ 
pretation follows from the simple circumstance that any 
source term in the energy balance equation of a general 
relativistic fluid may be formally rewritten in terms of an 
effective bulk viscosity [26} . 

The energy-momentum tensor of a general relativistic 
fluid with bulk viscosity as the only dissipative process 
can be written as 12611 


T^ v = {p + p + W)Uy.U v - (p + JTjg^, (57) 

where n is the bulk viscous pressure. The particle flow 
vector N M is defined as N M = nt/ M . In the framework of 


causal thermodynamics the entropy flow vector takes 
the form [38} 


= sN» - (58) 

where r is the relaxation time, and £ is the coefficient 
of bulk viscosity. In Eq. (l58l) , we have limited ourselves 
to considering only second-order deviations from equilib¬ 
rium. In the case of homogeneous and isotropic geome¬ 
tries, in the presence of bulk viscous dissipative phenom¬ 
ena, the energy conservation equation is obtained as 

p + 3{p + p + H)H = 0. (59) 

By comparing Eq. ES), giving the energy conserva¬ 
tion equation for a cosmological fluid in the presence of 
bulk viscosity, with Eq. (1341) , which includes in the energy 
balance the creation of particles from the gravitational 
field due to the curvature-matter coupling, it follows that 
these two equations are equivalent if 


Pc = n = - 


1 

3 H 


4t bv » ] 


2V' M - U'(4)L m 
U(4) 



(60) 

Therefore particle creation can be indeed described from 
a phenomenological point of view by introducing an effec¬ 
tive bulk viscous pressure in the energy-momentum ten¬ 
sor of the cosmological fluid. Hence it follows that the 
causal bulk viscous pressure n acts as a creation pressure. 

Hence, it would be interesting to investigate matter 
creation processes in modified gravity with a curvature- 
matter coupling from the point of view of bulk viscous 
thermodynamic processes. Hence we shall consider in the 
following that there is a change in the number of parti¬ 
cles, due to matter creation processes, with bulk viscous 
pressure playing the role of the creation pressure. We 
introduce a simple toy model in which the newly created 
particles obey, as a function of the particle number den¬ 
sity n, an equation of state of the form 


P = Po 



p={ 1 -l)p 1 (61) 


where po , no and 7 are constants, we have denoted k = 
Po/n 0 , and 1 < 7 < 2, respectively. Using Eq. (RJTl) . then 
Eq. (l59l) takes the form of a particle balance equation, 


n + 3Hn = Fro, 


where 


r = 


n 

7 H 


(62) 


(63) 


is the particle production rate, proportional to the 
bulk viscous pressure. Combining the equation of state 
Eq. (RJTl) with the Gibbs relation Tds = d{p/n) +pd(l/n) 
we obtain s = s " 0 = constant, that is, particles are created 
with constant entropy density. 





However, there is a major difference between the parti¬ 
cle creation irreversible processes in open thermodynamic 
systems and bulk viscous processes, and this difference 
is related to the expression for entropy production rate. 
While the entropy production rate associated to particle 
creation is given by 


= - 


3 Hp c 
T 



> 0 , 


(64) 


in the presence of bulk viscous dissipative processes the 
entropy production rate can be obtained as [26| 


V u ^ = - 


n 

T 


3 H 




n \3H- 


(65) 

In the particle creation model in open thermodynam¬ 
ics systems the entropy production rate is proportional to 
the creation pressure, while in the viscous dissipative pro¬ 
cesses thermodynamic interpretation is quadratic 

in the creation pressure, V fl S^ oc p 2 /£T, and, moreover, 
involves a new dynamical variable, the bulk viscosity co¬ 
efficient. 


IV. COSMOLOGICAL APPLICATIONS 


In the present Section, we consider several cosmologi¬ 
cal applications of the scalar-tensor formulation of mod¬ 
ified gravity with a linear curvature-matter coupling, as 
its interpretation as a particle creation theory. For a ho¬ 
mogeneous and isotropic geometry the gravitational field 
equations (071) - 11201) take the form 


3 H 2 = - + ( ^p+p 


3' 

v{ip) 


tp 


UW 

tp 


+ 2P'0 tp) - 2 U'{ip)L r , 


( 66 ) 


2H + 3H 2 = 


a ip 3 ip 


tp 


2P'(V>)-2 U'(iP)L n 


067) 


and 


- 3 (h + 2H 2 ) - V' {ip) + U' WO L m = 0, (68) 


respectively. 

Equation (1001) can be integrated to give the particle 
number time variation as 

n(t) = (69) 

cr 

where no is an arbitrary constant of integration, while, 
by assuming a barotropic equation of state of the form 
p = p (p) we obtain for the density evolution 


dp 


P + p{p) 


= —31na+ / T{t)dt + \npo, (70) 


where po is an integration constant. 


A. Cosmological models satisfying the condition 

217 W) - U'{ip)L m = 0. 


The linear curvature-matter coupling depends on the 
two arbitrary (and independent) potentials U{ip) and 
V{ip). As a simple toy model, in the following, we as¬ 
sume that the two potentials U and V are related to the 
matter Lagrangian via the relation 

2 V'{ip) - U'{ip)L m = 0. (71) 


Moreover, we restrict our analysis to the case of dust, 
with negligible thermodynamic pressure p = 0. With the 
choice of Eq. CD. the particle creation rate, given by 
Eq. (HIED, takes the form 

r=-kl,C Wi) ) = -kM^>0. (72) 

Therefore, in this approach the matter creation rate is 
determined by the function U{ip), which describes the 
coupling between the matter Lagrangian and the scalar 
field, its derivative, and the time variation of the scalar 
field ip only. 

With this choice, the variations of the particle number 
and of the matter energy density is given by 

^wv rti) = ? urn■ (73) 

With the assumption of Eq. CD on the potentials, the 
gravitational field equations (l66l) - (l68ll take the form 


^-V 2 r u -f-\ 


VW 

. tp 


- U'{ip)L n 


, (74) 


2 H + 3 H 2 


dip 1 U {ip) 1 V{ip) 

a ip 3^ ip 3 _ ip 

&(h + 2H 2 ^ = U'(ip)L m , 


U'{ip)L m , 

(75) 

(76) 


respectively. 

In the following, we consider only de Sitter type accel¬ 
erating solutions of the system given by Eqs. cam 
with H = H 0 = constant and apt) = exp {Hot). Then we 
obtain first 


U’ {iP) L m = 12 H 2 , 


(77) 


while the evolution equation for the scalar field ip, which 
can be obtained from Eqs. CD-CUb is given by 


ip + H 0 ip ■ 


Po 

with the general solution given by 


= 0 , 


(78) 


ip{t) = - 


1 


6H 2 


p 0 e~ 3Hot - 3 p 0 e~ Ho{t+2to) + 2 p 0 e~ 3Hot 


+6 H 0 iPoie- Ho{t - to) - 6H 0 {H 0 ipo + ip 0 i) 1 , (79) 
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where we have used the initial conditions ip (to) = ipo, 
and ip' ( t 0 ) = ip 0 i, respectively. Since for the dust fluid 
the matter Lagrangian is L m = p, Eq. C3 yields 


V W _ c 3Hot 


U(ip) 


P o 


or, equivalently, 


1 


dU(t) _ c 3H 0 t 


ip. 


(80) 


(81) 


U(t) dt po 

The above equation determines the time variation of 
the potential U(ip) as 

- 6 H 0 ip 01 e H °^ 2 _ 2Hn(t -^ 

Po 

(82) 

where Uq is an arbitrary constant of integration. The 
time variation of the particle creation rate is obtained as 

6H 0 U 0 


U ( t ) = Uo exp 


6 H 0 t 


m = 


Po 


2H 0 \ipoi\e H °( 2t+to) + p 0 e 2ffo( ‘- to) - Po 


(83) 

where we have assumed f/’oi < 0. The particle creation 
rate is non-negative for all times t, and is a monotonically 
increasing function of time for all t > to, and its initial 
value is given by T (to) = l2Ho\ipoi\/Po > 0. 

Therefore we obtain for the time variation of the mat¬ 
ter energy density the equation 


Po 

P(t) = jr exp 
Uo 


3e 2 H 0 (t-t 0 ) (2Ho\ipoi\e 3Hot ° + po) 


Po 


— 9H 0 t 
(84) ' 

For the time variation of the potential V ( t) we obtain 
the evolution equation 

dV dU 

~dt~ P ltt' ^ ^ 

which gives the potential V in an integral form as 

I" Hq \ 12 I V’Ol |e Jf ° (2t + t 0 ) _ 3 ^1 

V(t) = U 0 + 3H 0 U 0 el P0 \ x 

x [-p 0 e 2ffo(t - to) + 2Ho\ipoi\e H ° (2t+to) + p 0 ] dt, (86) 

where Vo is an arbitrary constant of integration. The 
creation pressure, defined as p c = —pT/(3H), is given by 

p c = 2U 0 [-poe 2 " 0 **-*") - 2ff 0 |^oi|e ffo(2t+to) + p 0 

’3e 2 H °(i-to) (2H 0 \ip 0 i\e 3H ° to + po) 


x exp 


Po 


-9 H 0 t 


(87) 


Finally, for the comoving entropy of the de Sitter type 
expanding Universe in the linear curvature-matter cou¬ 
pling theory we obtain 


S(t) = 


So 


So 

U (t) = U~o 6Xp 


+ 


3e 2H 0 (t-t 0 ) _ 6Hof 

6Hq \ipoi\ e H ° l 2t+t ° } " 


Po 


( 88 ) 


The entropy is a monotonically increasing function of 
time, with the property S(t) > 0, V£ > to- In the first 


order approximation, and for small times we obtain for 
the entropy the following expression 


S(t) oc exp 


/ 12g 0 2 |V>oi| 

V Po 


(89) 


Therefore, the curvature-matter coupling allows the 
production of a large amount of entropy during a de Sit¬ 
ter type evolutionary phase of the Universe. 


B. de Sitter type expansionary models with 
constant matter creation rate 


In the following, in the scalar-tensor representation of 
Lf (R, L m ) gravity, we consider a second simple cosmo¬ 
logical toy model by assuming that the cosmological ex¬ 
pansion of a dust Universe, with p = 0, is accelerating 
with a = exp (Hot), where Hq = constant, and the parti¬ 
cle creation rate is a constant during the entire accelerat¬ 
ing phase, and it is given by T = To = 3Ho = constant. 
Moreover, we take the matter Lagrangian as L m = p. 
Then from Eq. (1701) it follows immediately that the mat¬ 
ter density of the Universe is also a constant, 

P = Po = constant. (90) 


Then the constancy of the matter creation rate, given by 
Eq. (I5T1) imposes the following condition on the potentials 
U and V, 

2V + 3H 0 p 0 U = 0, (91) 

while the field equation Eq. (K?51l yields 

V — Po U = -6 H^Pp. (92) 

From the field equations Eqs. (1S51) and (1571) we obtain 
the evolution equation for ip as 

ip + H 0 ip + poU(ip(t)) = 0. (93) 

From Eqs. m and (1^71) we obtain 

^f^U + p 0 U = 6H 2 ijj. (94) 

By taking the time derivative of the above equation, 
and by eliminating ip with the help of Eq. (l93ll , it follows 
that U satisfies the equation 

U+^H o U+y h oU = 0, (95) 


with the general solution given by 


U(t) = 


e §H 0 (t a -t) 


215H 0 U 0 cos 


215-ffo i 

v / 215(5iL 0 U 0 + 4J7 0 i) sin 


-V215H 0 (t - to) 


- A VZH>Ho(t - to) 


,(96) 


where we have used the initial conditions U (to) = Uq 
and U' (t 0 ) = Uqi, respectively. The time dependence of 
the potential V is obtained from Eq. (15T1) in the form 
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V(t) = 


Poe 




860 Hq 


V215(2U 0 i - 19H 0 U 0 ) sin 


-V215H 0 {t -1 0 ) 


■ 43(5 H 0 U 0 + 2E/ 0 i) cos 


\V215H 0 (t -1 0 ) 


(97) 


With the use of the explicit time dependence of U in Eq. (1M1) . we obtain for ip(t) the following expression 
Ht) = 1161 1 Qg3 { - 430e H °(* 0-t ) (27^i - 3H oPo U 0 - 2 Po U 0 i) + 387 (30H^ o + 30H^ O i - 5H oPo U 0 - 2 Po U 01 ) 


+Poe i H o(to-t) |y 215 ( 39 H 0 Uo + UU 01 ) sin (^V215H 0 (t - t 0 )j + 43(15 H 0 U 0 - 2U 01 ) cos Q\/2l5 H 0 (t - t 0 ) 


(98) 


where we have used the initial conditions ip (to) = ipo and 
ip (to) = ipoi, respectively. 

Due to the cosmological particle production the en¬ 
tropy of the Universe increases as 

S(t) = S 0 e-T rWdt = S 0 e 3Hot . (99) 

However, the specific entropy s = S/V remains a con¬ 
stant during the cosmological evolution, s = so = 
constant. 


V. TOTAL ENTROPY BEHAVIOR IN Lf (R, L m ) 
GRAVITY WITH PARTICLE CREATION 

In the present paper, we have defined the entropy 
through the particle production rate, given by Eq. (43), 
as depending on the positive particle creation rate T via 
the relation 


In the standard thermodynamic description of physi¬ 
cal systems the time parameter t is not a thermodynamic 
equilibrium variable. Therefore, the variation of the ther¬ 
modynamic quantities should be considered with respect 
to some extensive variable. In the following we will adopt, 
following sun, as extensive variable for the cosmolog¬ 
ical system the proper volume enclosed by the apparent 
horizon, or, more specifically, its scale factor a. Then the 
relation between d/dt and d/da is simply 



( 101 ) 


In the following we denote by a prime the derivative with 
respect to the extensive variable a. 

We define the total entropy of a FRW Universe with 
dust as the sum of the entropy of the apparent horizon 
S a h, proportional to its area, and that of the matter par¬ 
ticles within it S m sun- For practical purposes, in the 
case of the flat FRW model, the total entropy is 


| = r > o. (ioo) 

Therefore, in an ever expanding Universe with particle 
creation, the matter entropy will increase indefinitely. On 
the other hand, all natural systems tend to approach a 
state of thermodynamic equilibrium, implying that the 
entropy of equilibrium systems never decreases, S > 0, 
and that it is concave when approaching the equilibrium 
state, S < 0. However, in the present, and several other, 
cosmological models, these fundamental requirements for 
the behavior of the entropy do not seem to be satisfied. 

The problem of the validity of the second law of ther¬ 
modynamics in cosmology was investigated in detail in 
la 123, where it was shown that the Universe approaches 
thermodynamic equilibrium in a de Sitter phase, if one 
defines the total entropy S tot of the Universe as the en¬ 
tropy of the apparent horizon plus that of matter and 
radiation inside it. Then it follows that Stot. increases, 
and that it is concave, thus leading to the result that the 
second law of thermodynamics is still valid for the case of 
the cosmological expansion. In the following, we investi¬ 
gate the thermodynamic properties of the total entropy 
in a Universe with matter creation. 


Stot — S a h + S m — + Qjjo n (t) , (102) 

where we have used the fact that the radius of the ap¬ 
parent horizon is r a h = H~ l {4(j. We also introduce 
an important observational quantity, the deceleration pa¬ 
rameter g, defined as 


_ d ! i _ H' i 


dt H 


H 


(103) 


Therefore, the variation of the total entropy can be 
obtained as 


J tot 


S'gh + S'm 

Sah + S rn 


(104) 


J tot 


q" 4 - q" 

J ah w O m 

Sah + S m 


(105) 


respectively. 

The particle number n satisfies Eq. ®, and is rewrit¬ 
ten as 


aH ——K 3 nH = Tn. 
da 


(106) 
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From its definition the derivatives of the entropy with 
respect to the scale factor can be evaluated as 


2t ra rr „_, , 47rn(a 

^tot uK ( a )“b 


H 2 { a y 


an = 2tt 

tot 3aH 4 (a) 


H 3 (a) [ 3 


F(o) 


— aH'(a) — H(a) 


(107) 


|2iJ(a) 3 a 2 H ,2 (a) 


+n(a) [a [T'{a) - 3 ciH"{a) + 12 H\a)} - 6r(a)] 
—3 H 2 (a) a [aH"(a) + H'(a)\ — 6 n(a) 

+2n(a) [T(a) - 3a£T »] 2 }, (108) 

respectively. 

In terms of the deceleration parameter we can express 
the variation of the total entropy with respect to a as 

4.7 t f I' \ 

S' tot = -(q + l) + 1 ^(qH+-)n, (109) 


S’L = (o)n(o) [T(o) - 3 aH'(a)\ 

+2H(a)n(a) a [r'(a) — 6q(a)H'(a)\ 3r (a)(q(a) — 1) 
—3 H 2 (a) a[q(a ) + 1 }H\a) + n(a) [ 6 g(a) — 2 aq'(a)\ 
+3a.ff(a)Y(a)}. (110) 


Therefore the standard thermodynamic requirements 
S'tot > 0 and S'tot — 0 impose the following constraints 
on the particle creation rate T, and its derivative with 
respect to the scale factor 


T(o)> 


3a [H(a) + 2n(a)] H'(a) 
2 n(a) 


+ 3 H(a), (111) 


T(a) > - 3[9{a) ^ ) H2{a) q(a)H(a), ( 112 ) 


and 


r'(a) < 


1 6 H(a) n{a) [a 2 H"(a) + 2T(a) — 4 aH\a)\ — a 2 H' 2 {a) — 2 n(a) [r(a) — 3aH'(a)\ 


2aH (a)n(a) 

+3 H 2 (a) [a ( aH"(a ) + — 6 n(a)] j 


(113) 


r (a) < 


— 2r(a)n(a) [T(a) — 3aH'(a )] + 3 H 2 (a) [a(g(a) + 1 )H'(a) + n(a) (6q(a) — 2a</(a))] 


2 aH(a)n(a) 

6 H(a)n(a) [T(a) + 2 aq(a)H'(a) + T(a)(— q(a))] — 3aff 3 (a)g , (a)|, 

I- 


+ 


(114) 


respectively. Since in the Lf ( R,L m ) gravity theory the giving 
matter creation rate T is determined by the coupling 
functions U and V, the thermodynamic conditions im¬ 
pose some strong constraints on the allowed physical form 
of these functions. 

A particularly interesting case is that of the de Sitter 
evolution of the Universe, with H = H$ = constant. For 
this situation the total entropy of the Universe is given 

by 


e 7r , 4?r 1 \ 

s “ = ni + m {) 


( 115 ) 


4ir / 47 r [F(a) — 3Ho] . . , . 

= -Tvvvn (a) = — 4 -- -n(a) > 0 , (116) 


tot 3 H$' 


3Hq 


S'toi = ^^{r 2 (a) + H 0 {aT\a) + 12H 0 \ 


-7H, 


.r(a)} 


< 0. 


(117) 


The thermodynamic condition of the non-negativity 
of the total entropy derivative with respect to the scale 
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factor imposes the conditions T > 3Hq and r'(a) < 
[7r(a) — T 2 (a)/H 0 — 12 iJ 0 ] /a on the particle creation 
rate T. In the particular case T = 3Hq, we obtain 
Stot = constant, showing that in this case the cosmolog¬ 
ical evolution is isentropic, with the total entropy being 
a constant. 


VI. DISCUSSION AND CONCLUSIONS 

In the present paper, we have considered the ther¬ 
modynamic interpretation of modified theories of grav¬ 
ity with a linear coupling between matter and geome¬ 
try, which we denote as Lf ( R , L m ) gravity. This theory 
represents a particular class, corresponding to a specific 
choice of the gravitational Lagrangian, of a very general 
class of theories, in which the action is an arbitrary func¬ 
tion of the Ricci scalar and of the matter Lagrangian. 
An interesting characteristic of these theories is the non¬ 
conservation of the energy-momentum tensor of the mat¬ 
ter, indicating that matter and energy fluxes can be gen¬ 
erated by the conversion of the geometric curvature, de¬ 
scribing the gravitational field, into matter. Hence the 
presence of matter, and its possible coupling to geom¬ 
etry could modify the cosmological evolution in a way 
that goes far beyond the standard description of gen¬ 
eral relativity. The presence of a source term in the 
energy balance equation can be naturally interpreted in 
the framework of the thermodynamics of open systems as 
describing a particle creation process, in which the “ge¬ 
ometric energy” of the gravitational field is transferred 
to “real” matter. During the particle production phase 
a large amount of entropy is produced. 

In order to estimate the effective thermodynamics 
quantities we have first introduced the equivalent scalar- 
tensor representation of the Lf ( R , L m ) theory, which can 
be formulated in terms of a scalar field ip, with two in¬ 
dependent potentials Viip) and U(ip), with the poten¬ 
tial U(ip) coupled to the matter Lagrangian. Using the 
scalar-tensor representation of the Lf ( R , L m ) theory, we 
have obtained the particle creation rate, the creation 
pressure and the entropy associated to the gravitational 
energy transfer to matter. The cosmological implications 
of the particle creation have also been investigated, by 
assuming a specific relation between the two potentials. 
The imposed condition makes the particle creation rate 
a function of the second scalar potential U(ip), which 
directly couples to the matter Lagrangian. The gravita¬ 
tional field equations corresponding to these choices have 
a de Sitter type accelerating solution, where the cosmic 
acceleration is triggered by the particle creation process, 
which generates a negative creation pressure. Thus, it 
was argued that the negative creation pressure is respon¬ 
sible for the accelerated expansion of the Universe. 

Matter creation processes are supposed to play a fun¬ 
damental role in the quantum field theoretical approaches 
to gravity, where they naturally appear. It is a standard 
result of quantum field theory in curved spacetimes that 


quanta of the minimally-coupled scalar field are created 
in the expanding Friedmann-Robertson-Walker universe 
[4lj |. That’s why finding an equivalent microscopic quan¬ 
tum description of the matter creation processes consid¬ 
ered in the present paper could shed some light on the 
physical mechanisms leading to particle generation via 
gravity and matter geometry coupling. In the following 
we will briefly point out that such mechanisms do ex¬ 
ist, and can be understood, at least qualitatively, in the 
framework of some semiclassical gravity models. 

In semiclassical gravity it is assumed that the gravita¬ 
tional field remains classical, while the classical bosonic 
fields f> are quantized. In order to couple quantized fields 
to classical gravitational fields the quantum energy mo¬ 
mentum tensor T^ v is replaced by its expectation value 
with respect to some quantum state T, thus leading to 
the effective semiclassical Einstein equation 14^ . 

~ 2 = c 4 (^1 1 ^) ■ ( H 8 ) 

Hence the classical energy-momentum tensor of the 
system is defined as |T) = T The semi¬ 

classical equation Eq. (111811 can be obtained from the 
variational principle |43[ 


6{S g + S^) = 0, (119) 

where S g = (1/167 tG) f Ry/^gd 4 x is the classical action 
of the gravitational field, and 


Sy = j Im^'i'll'^ - + a((^|^) - 1) dt, 

( 12 °) 

where H is the Hamiltonian operator of the system, and 
a is a Lagrange multiplier. The variation of Eq. (11191) 
provides the normalization condition for the wave func¬ 
tion ('k|'L) = 1, the Schodinger equation for the wave 
function 


( 121 ) 

as well as the semiclassical Einstein Eq. (11181) . In this 
simple case the Bianchi identities require the conserva¬ 
tion of the energy-momentum tensor, (T| T^ v |\I/) = 
0 . 

A very different set of semiclassical Einstein equations 
can be obtained by assuming a coupling between the 
quantum fields and the curvature of the space-time. In 
the model introduced in [43] the contribution to the to¬ 
tal action of the geometry-quantum matter coupling term 
was assumed to be of the form 


f RF((f(<P)))vV^d 4 x, ( 122 ) 


where F and / are arbitrary functions, and ((/(</>)))^, = 
('L(f)l f[(p{x)\ l'F(f)). Then, in the presence of such a 
geometry-matter coupling the Hamiltonian H(t) in the 
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Schrodinger Eq. (11211) is modified to [H} 


ritj), is governed by an equation of the form 


Hit) Hit) - INF ' «/($»* 

(123) 

where N is the lapse function, £* are intrinsic coordinates, 
such that the normal is everywhere time-like, and 7 = 
det "frs, where 7 rs is the metric induced on a surface 
a (f), which gives a global slicing of the space-time into 
space-like surfaces. The effective semiclassical Einstein 
equation takes the form jifjj 


h + 3 H (p+ + P4> ) + Fl {p * + P +> p * = 0, (125) 

where Ti is the particle decay rate, determined by the 
coupling between matter and geometry. For the energy 
density and pressure of the scalar field dark matter we 
can assume the standard form 

dp dp 1 

P<fi = Y + Uini (<t>), P<t,= Y~ U int(<t>), (126) 


Rfiv r, 




= 167tG 

— + g^OF 


(124) 


In Eq. (11241) the matter energy-momentum tensor is 
not conserved, ( 'T^ 0. Thus, this equation de¬ 
scribes an effective particle production process, and can 
be interpreted as giving an effective semiclassical descrip¬ 
tion of the quantum processes in a gravitational field. By 
modifying the classical part of the gravitational action we 
can recover the field equations Eqs. © and © , respec¬ 
tively, used in the present paper. Therefore the physical 
origin of the matter creation processes considered in the 
present paper can be traced back to the semiclassical ap¬ 
proximation of the quantum field theory in a Riemannian 
curved geometry. 

An interesting and important question is the physical 
nature of the particles that could be created via gravi¬ 
tationally induced creation processes. The most natural 
assumption would be that these particles are dark matter 
particles. It has been conjectured that dark matter may 
consist of ultra-light particles with masses of the order 
of m ~ 10 -24 eV (see (44j and references therein). From 
a physical point of view such a particle may represent a 
pseudo Nambu-Goldstone boson. Axions are other ultra¬ 
light dark matter candidates, with masses in the range 
m < I0 ~ 22 eV J4f|. Such extremely very low mass parti¬ 
cles can be created even in very weak gravitational fields. 
An alternative description of dark matter is provided by 
the so-called scalar field dark matter models [461]. in which 
it is assumed that dark matter is a real scalar field, mini¬ 
mally coupled to gravity, with the mass of the scalar field 
having a very small value of the order of m < 10 21 eV. 
For zero temperature scalar field dark matter models all 
particles in the system condense to the same quantum 
ground state, thus forming a Bose-Einstein condensate. 
Therefore scalar field dark matter models are equivalent 
to the Bose-Einstein condensate dark matter models [43] ■ 
This implies, from a physical point of view, that in the 
open irreversible thermodynamic model introduced in the 
present paper particle creation can take place also in the 
form of a scalar field. In such a model the evolution of 
the scalar field dark energy particles, with energy density 
p<f, and pressure p^, and having a particle number density 


where U((f>) is the scalar field self-interaction potential. 

The creation pressure corresponding to the scalar field 
creation processes can be obtained as 


(4>) _ Fi {Prfr + P<f>) Prp 

c SHn^ 


(127) 


It is interesting to note that Eq. (11251) . which describes 
the creation of a scalar field as a result of the geometry- 
matter coupling, can be written in an equivalent form 
as 


d> + 3Hd> + r (d>,d>,u^d> + u! nt {d)) = o, ( 128 ) 

where we have denoted T ^</>, d>, = Tip^/n^. There¬ 

fore in the scalar field dark matter model a friction term 
in the scalar field evolution equation Eq. (11281) does ap¬ 
pear naturally, and in a general form, as a direct conse¬ 
quence of the irreversible thermodynamics of open sys¬ 
tems as applied to the dark matter case. Hence scalar 
field dark matter can be a result of the cosmological par¬ 
ticle production due to the geometry-matter coupling in 
modified gravity theories. For gravitational models with 
an action given by an arbitrary function of the Ricci 
scalar, the matter Lagrangian density, a scalar field and a 
kinetic term constructed from the gradients of the scalar 
field, respectively, see j4S|. 

The Lf(R,L m ) gravitational theory investigated in 
the present paper predicts the possibility that matter 
creation, associated with the curvature-matter coupling, 
could also occur in the present-day universe, as proposed 
by Dirac [HI a long time ago. The late expansion of 
the Universe 0,1 may be considered as an empirical 
evidence for matter creation, and a viable alternative 
to the mysterious dark energy. Presently the existence 
of some forms of the curvature-matter coupling lead¬ 
ing to matter creation processes cannot be fully ruled 
out by the existing cosmological observations or by as- 
trophysical data. Presumably, the functional forms of 
the potentials U(i/>) and U{i/j) that completely charac¬ 
terize the Lf (i?, L m ) gravitational theory will be pro¬ 
vided by fundamental quantum field theoretical models 
of the gravitational interaction, thus opening the possi¬ 
bility of an in depth comparison of the predictions of the 
Lf (R. L m ) gravity with cosmological and astrophysical 
observational data. 
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